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1. INTRODUCTION 

The Electro Cardiogram (ECG) produces electrical signals of the each cardiac cycle. In the cardiac 
cycle,each event has its own significance to study the behaviour of patient cardiac pathophysiology, 
Generally ECG signals are Bio electrical signal which gives the electrical activity of heart versus time. 
Therefore it is very important to diagnose for analysing heart function [1]. The Bio electrodes are placed on 
the skin of the patient to acquire ECG signals. The pacemakers are located in the upper part of the right 
atrium. It fires electrical pulses to the nerves to stimulate the contraction phase. These pulses extend over the 
atrial walls and activate cardiac muscles to contract. 

The ECG frequency range is typically 0.05 to 100Hz and power-line signal’s frequency range is 
50Hz. So, ECG signals sensitive to the power line signals in the range around 50Hz which are causing 
interference [2]. 50Hz PLIN will interrupt the P and Q waves of the ECG signal. Most of Asian regions, 
Domestic and hospital power line are in the range of 50 Hz. So the frequency components associated in ECG 
that is 50 Hz are effected by the power line signals causing interference in ECG. But lack of power line 
quality the power signals are swings between the 47 to 53 Hz so the interference also effects from this range 
of power line signal. To mitigate this dynamic interference from power line we need use adaptive filter to 
suppress this random noise causing from power line [3]. 

An adaptive noise elimination filter has been used to evade this impending loss of information. Four 
different waves can be observed while recording ECG signal those are PQRST. The depolarization of right 
atria represents P wave. While the rapid depolarization of right and left ventricles represents QRS wave. The 
repolarization of the ventricles represents T wave [4], [5]. Any deviation in the said parameters leads 
abnormalities in the heart .The wave form related to QRS complex represent the contraction of left and right 
ventricles, which is more powerful than that of atria .It comprises muscle mass and causing a more ECG 
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deflection. The Q wave signifies the signal horizontal (i.e. left to right) current as a potential travel through 
the inter-ventricular system. The Q wave is not having a septal origin shows myocardial violation which 
involves the full depth of myocardium. 

The P wave arises when the SA node (Sinus Atria) generates a potential which depolarizes the atria. 
However as long as the atrial depolarization takes place to spread through the AV node to the ventricles, each 
P wave should be trailed by QRS complex [6]. From the commencement of QRS complex is called PR 
interval. This indicates the time that it takes for the electrical impulses produced in the SA node and to travel 
through the atria and across the AV node (Atria Ventricle) node to the ventricle.In adaptive filter, least mean 
square algorithm requires input signal and reference signal to update the (tap weights) filter coefficients of 
the adaptive FIR filter. For every iteration, LMS algorithm update new tap weights based on the previous tap 
weights to minimize the error. After several iterations, it eliminate the noise in the adaptive filter and gives 
best minimum mean square error. This method follows the computation based on the past available 
information. The RLS (Recursive least square) algorithm gives better convergence than LMS 
algorithm [7]. The main disadvantage of the RLS algorithm is having high computational cost. In the paper, 
section II shows the LMS and RLS algorithm, section III gives the simulation results and section IV gives the 
conclusion results. 


2. RESEARCH ADAPTIVE FILTERING 

Adaptive filter involves the alteration of filter parameters (coefficient) over time to reduce the noise 
in the signal and to minimize the error [8]. Digital signal processing exhibited by most of the adaptive filters 
will be digital in nature because of complexity in optimizing algorithms. Adaptive filters are best suited when 
there is large uncertainty and filter has to compensate that or signal conditions are slowly changing. The 
performance of adaptive filter involves two process,which are filter processing and adaptation process [9]. 
The adaptive filter output Z(n) is given as (1). 


Z(n)=P(n)S(n) (1) 


Where S(n) is the input ECG signal and P(n) are the adaptive filter coefficients. 

As we known ECG signal interference by 50 Hz power line signal .so, we have to suppress the 50 
Hz component in the ECG signal which are originated by the power line. If we remove complete 50 Hz 
component in the ECG signal there may be data loss in the ECG signal which are associated with 50Hz 
frequencies region. So we need to estimate the disturbance by power line signal for that we giving power line 
signal as reference to the adaptive filter. Now the adaptive filter will track the power line components in ECG 
signal and we can easily extract interference part in the ECG signal. The difference between the desired 
signal d(n) and the signal from the output of the filter Z(n) is the error signal e(n). 


e(n) = d(n) — Z(n) (2) 
The filter variable updates filter coefficients at every time instantaneous. 
P(n + 1) = P(n) + AP(n) (3) 


Where P(n+1) is the updated weight vector with the previous weight vector and A is the correction factor 
depends on the value. 


2.1. Least mean square algorithm 

These algorithms was suggested by Widrow and Hoff. They developed LMS from their studies of 
pattern recognition [10]. In the LMS algorithms, the correction applied to the above mentioned estimate 
includes product of three factors: the error signal e (n-1), the (scalar) step-size parameter (u) and the 
tap-input vector s(n-1).LMS algorithm is the best selection if we are dealing with adaptive digital circuits in 
this case a filter that reject the bands of signal those causes interference to the ECG [11]. The basic operation 
of LMS algorithm is the recursive updating nature of filter coefficients by the reference of error signal. Each 
iteration of LMS involves three steps: 
Filter output: 


Z[n] = Ynzo S[n]P[n] (4) 


Estimation error: 


Int J Elec & Comp Eng, Vol. 9, No. 1, February 2019 : 505 - 511 


Int J Elec & Comp Eng ISSN: 2088-8708 Oo 507 


e(n) = d(n) — Z(n) (5) 
Tap-weight adaptation: 
P[n + 1] = P[n] + pS[nJe[n] (6) 


d(n) is taken as desired signal and Z(n) is the output response of the adaptive filter equation 4 shows the 
output response of the Filter with input signal S(n) and filter coefficients P(n). 


2.2. Recursive least square algorithm 

With supreme computational complexity, RLS is the fastest converging algorithm. It cancels 
maximum amount of noise by minimizing error with the fastest rate. So in this study, a tradeoff between 
computational complexity and convergence rate is done to attain the utmost noise free signal. The filter 
output and error function of RLS algorithm are shown in (8) and (9). 


ATIL(n-1)u(n-1) 
1+A71 uHL(n-1)u(n) 


R(n) = (7) 


Where R(n) is the vector gain, L(n) is the inverse correlation matrix, u(n) is the buffered input vector and 
A7* denotes the reciprocal of the exponential weighting factor. The filter output is 


Z(n ) = P™(n — 1)u(n) (8) 


Error signal: 


e(n) = d(n) — Z(n) (9) 
The updated coefficients as shown in (10-12) 
P(n) = P(n — 1) + R” (n)e(n) (10) 
P(n) = P(n — 1) + R" (n)[d(n) — Z(n)] (11) 
P(n) = P(n — 1) + R" (n)[d(n) — PT(n — 1)S(n)] (12) 


where R" (n) is the gain constant, with a sequence of training data up to time, the RLS algorithm 
estimates the weight by minimizing the resulting cost [12]. Where u(n) is the input, and A is the stabilization 
parameter. 


3. RESULTS AND ANALYSIS 

ECG signals are randomly taken from MIT BIH data base as shown in Figure 1 i.e. 
(101,104,106,109,124) [13-15]. The length of each ECG signal is restricted to 3600 samples. In LMS 
algorithm the mu is taken as 0.02 [8] and in RLS algorithm lambda is taken as one randomly. Matlab is taken 
as a tool for simulations and noise is considered as Adaptive White Gaussian Noise (AWGN) with power line 
interference of 50Hz on ECG signal with less noise power. After de-noising using LMS and RLS algorithm, 
ECG signals peaks are estimated and results are compared with original signals without noise. 

Implementation of LMS Algorithm by reducing channel Noise and reducing Power line interference 
of 50 Hz in ECG Signal. The Figure 1 shows the ECG signal 101 from the database, Figure 2 is noisy signal 
i.e AWGN. Figure 3. is corrupted ECE signal from noise. Figure 4. indicates the processed ECG signal after 
suppression of AWGN using LMS algorithm with Mu is taken as 0.02. Figure 5 specifies Power Line 
Interference suppression using LMS Algorithm with Mu is taken as 0.02. Figure 6 indicates the processed 
ECG signal after suppression of AWGN using RLS algorithm with lamda is taken as one. Figure 7 shows 
Power Line Interference suppression using RLS Algorithm.Different patient signals 101, 104, 106, 109, 124. 
Wav files are taken from MIT BIH database. The signal power and noise power of ECG signals are measured 
before applying to the adaptive algorithms. 
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Figure 3. ECG signal+noise Figure 4. Noise removal in ECG signal 
using LMS algorithm 
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Figure 5. (50Hz) Power line interference suppression Figure 6. Noise removal in ECG signal 
using LMS algorithm using RLS algorithm 
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Figure 7. (50Hz) Power line interference suppression using RLS Algorithm 


The convergenge rate is not same for LMS and RLS algorithm [8]. These signals are processed 
through Severval iterations by LMS and RLS algortihms in the adaptive filters to minimize the error. Table 1 
and Table 3 estimates the power of the error signals and mean square deviation. The same signals are taken 
with 50Hz power line interference by considering uniform noise on each signal. It is observed that more the 
signal power, minimizes the error as shown in Table 2 and Table 4. Table 5 shows that ECG signals under 
noiseless conditions, QRST peaks are estimated. Table 6 and Table 7 shows that ECG signals corrupted with 
noise and powerline interference, QRST are measured after processed through adaptive filter using LMS 
algorithm. It is observed that QRST peaks are influenced with noise, but the heart beat rate information is 
appropriate.The same case is taken using RLS algorithm to measure QRST peaks with more influenced with 
noise. It is observed that heart beat rate information is appropriate and removals on QRST peaks are 
minimized. By observing Table 5 to Table 9 power line interference (50Hz) in the system and noise in the 
channel are suppressed to maximum extend using LMS and RLS algorithm. The peaks of QRST in ECG 
waves and heart rates are measured for different signals from MIT_BIH data base. 


Table 1. Channel Noise Reduction on ECG Signal using LMS Algorithm by Considering Interference of the 
Same Signal by ith Different Phase in the Channel with Mu=0.02 Taken Randomly 


Signals from ECG Signal Noise Signal Power of the Error Mean Square 
MIT BIH data base Power (dB) Power (dB) Signal (dB) deviation (dB) 
101 -0.3339 0.0843 -8.9620 -21.2726 
104 -0.1954 0.1034 -8.7943 -21.0287 
106 -0.0922 0.0702 -8.6967 -21.0243 
109 -0.4866 0.0771 -9.1009 -21.0891 
124 -1.2148 -0.032 -9.8297 -22.5220 


Table 2. Power Line Interference of 50Hz on ECG Signal by Reducing Interference using LMS Algorithm by 
Considering Interference of the Same Signal by the Different Phase in the Channel 
Types of Signal from ECG Signal Noise Signal Power of the Signal Power of the Error 


MIT BIH data base Power (dB) Power (dB) and Noise (dB) Signal (dB) 
101 -0.3339 -3.010 5.6820 -2.9158 
104 -0.1954 -3.010 5.8206 -2.7737 
106 -0.0922 -3.010 5.9218 -2.6680 
109 -0.4866 -3.010 5.5267 -3.0627 
124 -1.2148 -3.010 4.8001 -3.7882 


Table 3. Channel Noise Reduction on ECG Signal using RLS Algorithm by Considering Interference of the 
Same Signal by ith Different Phase in the Channel with A =1 Taken Randomly 


Signals from MIT-BIH ECG signal Noise signal Power of the 
Arrhythmia data base Power (dB) Power (dB) errorSignal(dB) 


101 -0.3339 -0.0467 -0.0274 
104 -0.1954 -0.0472 -0.0581 
106 -0.0922 -0.1854 -0.0630 
109 0.4866 -0.2264 0.0268 
124 -1.2148 -0.0494 0.0065 
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Table 4. Power Line Interference of 50Hz on ECG Signal by Reducing Interference using RLS Algorithm by 
Considering Interference of the Same Signal by the Different Phase in the Channel 


Signal from MIT- 


Power of the original 


Noise signal 


Power of the Signal 


Power of the 


BIH data base ECG signal (dB) Power (dB) and Noise (dB) Error Signal (dB) 
101 -0.3339 -3.0103 5.6879 -0.0274 
104 -0.1954 -3.0103 5.8274 -0.0581 
106 -0.0922 -3.0103 5.9320 -0.0630 
109 -0.4866 -3.0103 5.5360 0.0268 
124 -1.2148 -3.0103 4.8114 0.0065 


Table 5. Each ECG Signal’s Length is Having 3600 Samples Under Noise Less Conditions, QRST Peaks are 
Estimated in ECG Signals as Shown in the Table (mV) 


E ane A R Peak (mV) SPeak (mV) TPeak (mV) Q Peak (mV) Re cae 
101 203.23 "25.535 35,965 6.672 45.83 
104 266.12 3.458 259.47 141.616 58.33 
106 355.24 -26.902 145.803 115.92 54.16 
109 256.372 23.832 208.24 163.611 66.66 
124 390.59 -41,245 20,3103 102.8 33.33 


Table 6. LMS Algorithm is Taken on Noisy ECG Signals which are Corrupted with Additive 
White Gaussian Noise. QRST Peaks are Measured (mV) 


re ene R Peak (mV) $Peak(mV) TPeak (mV) Q Peak (mV) oe ea 
101 75.65 -9.469 13.200 -2.706 45.83 
104 98.99 1.317 96.503 52.671 58.33 
106 132.36 -9.690 54.309 42.87 54.16 
109 94.71 8.86 77.436 60.49 66.66 
124 145.15 -15.42 7.398 38.233 33.33 


Table 7. LMS Algorithm is Introduced on Signal of ECG to Minimize the 50Hz Power Line Interference on 
ECG System for Considering mu=0.02.The Readings are Taken in Terms of Milli Volts 
as Shown in the Table 


pesca MIBI R Peak(mV) S$ Peak(mV) © TPeak(mV) QPeak(my) {Heart Rate 
101 75.745 9.571 13.242 2392 45.83 
104 98.855 1.691 96.343 52,503 58.33 
106 132.288 -9.827 54.395 42.822 54.16 
109 96.264 8.813 77.511 61.396 66.66 
124 144.994 -15.207 7.584 38.137 33.33 


Table 8. RLS Algorithm is Taken on Noisy ECG Signals to Reduce AWGN Noise on the Noisy ECG Signals 
Readings are Taken for Different Signals inMIT Data Base and Measurements are in Milli Volts 


etna He ais R Peak (mV) SPeak(mV) TPeak (mV) Q Peak (mV) eei eg 
101 0.208 -0.027 0.036 -0.003 45.83 
104 0.266 0.003 0.270 0.140 58.33 
106 0.341 -0.029 0.144 0.112 54.16 
109 0.261 0.027 0.214 0.161 66.66 
124 0.446 -0.047 0.024 0.115 33.33 


Table 9. RLS Algorithm is Considered to Remove Power Line Interference on the ECG Signals Which are 
Interfered with 50Hz Power Line Frequency in the System Volts 


Sota ts — R Peak (mV) SPeak (mV) TPeak (mV) Q Peak (mV) ae fol 
101 0.2089 -0.0275 0.0371 -0.0035 45.83 
104 0.2660 0.0032 0.2702 0.1406 58.33 
106 0.3412 -0.0299 0.1447 0.1056 54.16 
109 0.2623 0.0272 0.2142 0.1611 66.66 
124 0.4459 -0.0472 0.0246 0.1154 33.33 
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4. CONCLUSION 

In this paper, the ECG signals are taken from physio-net for analysis. For Bio-telemetry 
applications, signals are transmitted from remote locations, LMS and RLS adaptive algorithms are 
considered for suppression of AWGN noise and power line interference on ECG signals. The analysis is 
carried out to evaluate the ECG signal power, noise signal and Mean square error. It is observed that Mean 
square error is less in RLS then LMS. The rate of convergence is more fast that LMS algorithm. Noise 
cancellation capacity is good in RLS than LMS. But the implementation is complex over LMS.After 
denoising at the receiver end, QRST peaks and heart beats are estimated and compared with the original 
QRST peaks and heart beats. 
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